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Abstract 

Spontaneous directed motion, a hallmark of cell biology, is unusual in classical statistical physics. 
Here we study, using both numerical and analytical methods, organized motion in models of the 
cytoskeleton in which constituents are driven by energy-consuming motors. While systems driven 
by small-step motors are described by an effective temperature and are thus quiescent, at higher 
order in step size, both homogeneous and inhomogeneous, flowing and oscillating behavior emerges. 
Motors that respond with a negative susceptibility to imposed forces lead to an apparent negative- 
temperature system in which beautiful structures form resembling the asters seen in cell division. 
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Spontaneous directed motion driven by active processes is crucial to biology. Such mo- 
tion is only possible because the cell is a far-from-equilibrium many-body system. The 
cytoskeleton of eukaryotic cells is built, maintained and adaptively reorganized through ac- 
tive transport and force generation powered by ATP hydrolysis. Oscillations of the mitotic 
spindle during cell division^ and cytoplasmic streaming^ dramatically illustrate that the cell 
is not at equilibrium. Driven motions of cells are also important at higher levels of orga- 
nization in living things ranging from mechanosensation^ to the developmental processes in 
which the genetic code unfolds to create a multicellar organism^. Sustained spontaneous 
collective motion is quite remarkable in many-body physics. Superfluidity and supercon- 
ductivity are examples of metastable states of motion made possible by quantum statistics. 
The biological example provided by the cytoskeleton is seemingly quite different leading 
not to infinitely long-lived states but to ones that go away when the cell is depleted of fuel 
and dies. Nevertheless, like the quantum examples, the motion of the cytoskeleton is an 
emergent many-body phenomenon reflecting broken symmetries. 

Here we explore the origin of spontaneous collective motion for systems of many interact- 
ing biomacromolecules with motor-driven active processes using a systematic perturbative 
expansion of the many-body master equation treating nonequilibrium motorized processes. 
We model the motors as generating a time series of isotropic kicks on the constituents of 
a many-body assembly. Earlier- we showed that quite generally the corresponding master 
equation, when expanded to the lowest order in the kick step size, yields an effective tem- 
perature, T e ff, which explicitly depends on the total motor activity and on the way in which 
motors respond to imposed forces. A system described by an effective temperature alone^ - - 
cannot undergo spontaneous directed motion unless it is quantum mechanical. Pursuing the 
expansion to higher order, however, reveals the possible emergence of spontaneous directed 
collective motion quite generally from a quiescent homogeneous state, albeit one with rigidity 
owing to broken translational symmetry, as in a glass. The underlying dynamic instability 
is induced by a sufficiently strong internal agitation in terms of kick step size. This provides 
a general mechanism for spontaneous flows in an active assembly of interacting constituents. 

Combining a linear stability analysis with a trial solution of the many-body master equa- 
tion allows us to identify possible dynamic phases that depend on the motor kick step size 
and susceptibility. We find that for sufficiently large kicks and high activity, susceptible 
motors, i.e., motors whose kick rate depends on the forces exerted on them, can generate 
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spontaneous flow, whereas adamant motors, indifferent to imposed forces, would merely 
drive fluidization of an active system. We have also carried out simulations on a minimal 
cytoskeleton model incorporating motor dynamics to compare with our analytical predic- 
tions. The simulations not only verify the predicted phase diagram, but also highlight how 
the combination of network connectivity with motor susceptibility determines the formation 
of nonequilibrium structures. The simulations show an oscillatory phase separation at inter- 
mediate network connectivity and formation of aster-like patterns/bundle-connected poles 
when driven by motors with negative susceptibility, i.e., motors that move against the force, 
energetically uphill. The latter corresponds to a negative-temperature system where inter- 
esting structures emerge much like vortex condensation in two-dimensional turbulence*^. 

We are far from the first to try to understand the physics of spontaneous collective motion 
in biology. Jiilicher and Prosit studied a one-dimensional stochastic model which assumed 
an underlying ratchet potential already breaking translational symmetry. Motor coopera- 
tivity then leads to a dynamical phase transition to spontaneous directed motion despite the 
system's spatial symmetry. Thinking of the cytoskeleton, an assembly of filamentous polar 
polymers actively connected by crosslinkers, as an active polar gel has allowed the construc- 
tion of continuum theories, based on conservation laws and symmetry considerations, which 
also generate active flows^— . Pattern formation in active fluids has also been discussed 
based on a reaction-diffusion-advection mechanism^. 

Here we model the stochastic nature of the motor kicking via a master equation for the 
many-body probability distribution function \l/({r},t)^^ 

J^({rM) = (L FP + L NE )M/({f},t). (1) 

Here L FP = D £\ V.; ■ Vj — D /3 J2i ' ( — Vji7) is the usual many-body Fokker-Planck oper- 
ator describing passive Brownian motion with Dq denoting the ordinary diffusion coefficient 
at ambient temperature T and /3 = l/k^T. The gradients of the many-body interaction 
free energy U({r}) = U(ri,r 2 , ■ ■ ■ ,r n ) = J2<ij> u (^i ~~ ^j) & Ye the local forces acting on 
individual particles, where r*j is the position of the ith particle and (• • • ) denotes the nearest 
neighbor pairs. The effects due to nonequilibrium motorized processes are summarized by an 
integral kernel L NE ^({f},t) = / lUdP^K ({?} -»■ {r})V({P},t) - K({f} -»■ {P})^({r}, t)], 
where K({r'} — > {r}) encodes the probability of transitions between different particle con- 
figurations per unit time. Motor kicking noise is a finite jump process with a rate that 
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depends on whether the free energy is increased or decreased when a step is made 

k = k[Q(AU) exp(-s u /3AU) + Q(-AU) exp(-s d /3AU)}. (2) 

Here G is the Heaviside step function and AU = U + l \ —U (r) is the energy change due to 
the kick identified by a vector I = In. The kick step size I and the basal kicking rate k define 
the dimensionless motor activity A := k,1 2 /D , an analog of the Peclet number in turbulent 
diffusion. This model rate couples the chemical reactions leading to the motor activity to the 
local mechanical forces acting on the motor being parametrized by the susceptibility s which 
may take different values for uphill (s u ) moves and for downhill (sd) moves depending on 
the biochemical mechanism of the motors. When s — > 1 the motors are susceptible, slowing 
down when they climb up against obstacles and accelerating when they move energetically 
downhill; in contrast s — > corresponds to completely adamant motors which kick at a rate 
unperturbed by the free energy landscape. 

To examine the small kick limit, we first expand the equation in powers of I up to the 
quadratic order. The simplest case, isotropic kicking and symmetric susceptibility (i.e., 
s u — Sd = s) leads directly to an effective Fokker-Planck equation^ 

^m,t) = D eS J2{^ - V, • [(-ViArfrZT)*] }, 

i 

where 

/ 1 Kl 2 \ 

D °« = D °{ 1 + 2dD-J> 

These simple expressions ©-(E]) valid for general spatial dimensions d have nontrivial im- 
plications. In the small kick limit, the active system, while out of equilibrium, behaves as if 
it is at an effective canonical equilibrium characterized by an effective temperature T c g. The 
effective diffusion constant D e ^ (equation 2]) is enhanced by the active processes regardless 
of motor adamancy, consistent with recent observations of enhanced cytoplasmic diffusion^. 
T c g (equation [5]) is fully determined by the motor activity A = k1 2 /Dq and the motor suscep- 
tibility s; susceptible motors with s > 1/2 yield T e g < T. When motor activity dominates 
over thermal noise, i.e., A ^> 1, the effective temperature diverges as T e g/T ~ l/(2s) as 
s — > 0. Thus intense kicking by adamant motors leads to a very high effective temperature 
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just as observed in experiments^ and simulation studies 2 ^ 2 ^. A more detailed discussion 
can be found in a separate work-. 

To probe the dynamic instability that may give rise to the spontaneous motion, we must 
go beyond the effective equilibrium and expand to quartic order of the kick step size obtaining 

^m,t) = D eS j2{^-^i-[(-v^ eS um} 

i 

+ Kl A (cos A 9) ri x ^F 4 (vl m) f/, vS n) ^). (6) 

i 

The complicated functional Fi is the divergence of a flux, i.e., Fi = — Vj ■ J", where Jf is 
the probability current due to active events on particle % given by 

- J? = + ^V, (V 2 /3C/*) + |v,/H7V?¥ 

2 3 

+ S -(V t pU) 2 V^ + 8 -{V t puf^- (7) 
4 o 

Therefore while at quadratic order in I a motor-driven system exhibits enhanced diffusive 
dynamics at an effective equilibrium, at quartic order, a net streaming flow becomes possible. 

In an early study of the stability and dynamics of a motorized assembly Shen and 
Wolynes^ pictured the motors as introducing a modification to the Debye- Waller factors 
of the localized particles. They found an expression of the deviation of the total local- 
ization strength a from its thermal value a in terms of the motor properties. Thermal 
self-consistent phonon theory— gives a for a central particle that depends on the a of all its 
neighbors. Combining these two aspects allows a self-consistent determination of mean-field 
(a, a) solutions allowing an identification of stability limits. 

Assuming s u = Sd = s, the second moment closure^ reduces to a simple expression 
(a — a) I a — (s — 1/2) exp [s(s — l)al 2 } K,l 2 /dD . Thus for s = 1/2 (T e ff = T) chemical 
noise does not modify the mechanical stability (a = a); for s < 1/2 (T e g > T) stability is 
weakened (a < a) whereas for s > 1/2 (T c g < T) stability is enhanced (a > a). 

For spontaneous collective motion to occur there must be a nontrivial dynamic first 
moment that indicates a moving fiducial lattice. The second moment still has its steady- 
state value describing vibrations about the fiducial configuration. We thus write down a 
trial function of the master equation as a collection of Gaussians with moving centers and a 
steady variance 

({Fi}; d) = IT (a/ir) d/2 e -^-^(t)]\ ( 8 ) 
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where Wi — fi — Ri denotes the displacement of particle % from its equilibrium position 
fhi(t) = (wi\ty) defines the dynamic first moment of particle i, and the total localization 
strength a of individual particles is inversely related to the second moment. Applying the 
first moment closure d t {wi\^>) = (wi\(L F p + L NE )\1/) to equation dSJ) with this ansatz (jSJ) as 
well as the coupled-oscillator expansion of the effective potential leads to coupled equations 
for the mi's. 

To investigate the spontaneous emergence of directed motion from a quiescent homoge- 
neous state, we carry out a linear stability analysis on the non-moving state (i.e. rhi — 0). 
In view of the biological relevance of collective motion in one-dimensional scenarios, such 
as the filament sliding in motility assay^ and flow in the cell cortex—, we first demonstrate 
the instability onset for the ID case. Detailed analysis for general dimensions will be given 
later. Consider a spatially varying trial solution of the form 

m i = rhe rkt+ikR \ (9) 

where m denotes the amplitude of the first moment, k the wavenumber of the spatial modula- 
tion of a collective mode and the growth rate of the fc-mode. Even if kicks are isotropic and 
the interaction Hamiltonian preserves rotational symmetry, spontaneous symmetry breaking 
occurs giving flow in a specific direction. 

The linearized equations for the first moments in the long-wavelength limit (expanded to 
quadratic order in k) are 

dtmi = y|^^-)(Jij - i^) 2 W S ;a,a)mi, (10) 
j 

where 

f(l,s;a,a) = —Dq — snl 2 + s 2 Kl 4 a 

Since dtmi = r^mi, the growth rate is proportional to k 2 f(l, s; a, a) up to 0(k 2 ). It follows 
that a strictly uniform state (k = 0) would not undergo small-amplitude dynamic instability 
regardless of the motor activity and susceptibility but the nonmoving state is barely stable in 
the absence of spatial modulation. Moreover, the sign of /(/, s; a, a) determines the stability 
behavior for finite k. When the kick step size I is small, / is negative indicating diffusive 
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relaxation towards the non-moving state. As / increases, i.e., as the dynamics is more 
dominated by active processes, instabilities grow: /, thus r^, becomes positive, signifying 
spontaneous initiation of a collective flow in a spatially modulated state, when I exceeds a 
threshold value Z t h given by 



Consider the limit of high motor activity, such that a ^> a and k ^> aD^, we then have 
~ l/(sa) + Dq/(sk), indicating that high motor susceptibility s and kicking rate k as 
well as large mechanical gradient a lead to a low instability threshold in terms of kick 
step size (see Supplementary Fig. SI for a detailed illustration). For the case of asymmetric 
susceptibility (s u ^ sj), the factor (s u — s^) accompanies all the cubic-and-above odd powers 
in I in the expansion, leading to a lower threshold for the kick size compared to that for the 
symmetric case where deviations from effective equilibrium start at quartic order in /. 

We carried out the self-consistent calculation described earlier on a minimal model of 
the cytoskeleton as a cat's cradle^ to determine (a, a). We now use this to obtain the 
growth rate. As before^, we study a model network consisting of nonlinear elastic filaments 
characterized by relaxed length L e and stretching stiffness /3j built on a three dimensional 
random lattice of crosslinks at density p. The network connectivity P c is defined as the 
fraction of nearest-neighbor pairs of crosslinks connected by filaments. The unit of length is 
the average separation between the neighboring crosslinks. 

In Fig. Q] we show the phase diagrams for possible dynamic states as a function of kick 
size / and susceptibility s for several values of the network connectivity. In all the cases, 
there are two stability boundaries, one for small s (s < 1/2), one for large s (s > 1/2). In 
the low-s regime, as I reaches a critical value, l*{s) (blue dotted line), finite solutions for 
(a, a) become unstable, the vanishing of the localization strength then suggests the system 
becomes fluidized. In the high-s corner, as the kick step size exceeds a threshold value 
/ t h (lower boundary of the color-coded region), an exponential growth of the first moments 
occurs for small but finite k modes. This indicates the emergence of modulated flowing 
states. In this region, stable finite (a, a) solutions exist with a being considerably larger 
than a, reflecting the enhancement of stability by susceptible motor kicking. Note that 
as motor susceptibility increases, the threshold kick step size decreases. In the rest of the 
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FIG. 1: Phase diagram for possible dynamic states. The parameter plane indicates the motor 
kick step size I and the motor susceptibility s. The model cytoskeleton used to obtain a and a is 
characterized by crosslink density p = 0.8, relaxed length of the filaments L e = 1.2 and stretching 
stiffness /3-y = 5. From panel (a) to (c), the network connectivity varies with P c = 0.2,0.5, and 
1, respectively. Dq = 0.1 and k = 20. In the flowing regime, there are stable nontrivial (a, a) 
solutions and positive r^; in the diffusive/relaxational regime, while there are stable nontrivial 
(a,6t) solutions there is a small negative r^; in the fluidized regime, finite (a,a) solutions are 
unstable. As network connectivity rises, the flowing phase region expands whereas the fluidized 
state region shrinks. The logarithm of the normalized growth rate rk/k 2 for flowing instability is 
color-coded, showing the increase of instability with I and s. 



diagram, a and a are comparable and the negative growth rate indicates diffusive modes. 
(Close to detailed balance, s = 1/2, diffusive modes persist over the entire relevant range of 
I) 

In the figures, we color-code the logarithm of the normalized growth rate rk/k 2 for the 
flow instability; the growth rate increases with kick size I and susceptibility s. Comparing 
the diagrams for different values of connectivity, we see that as P c increases, the region 
corresponding to the fluidized state shrinks, since increasing the number of bond constraints 
stabilizes the system against fluidization. On the other hand, the region corresponding 
to flow expands toward lower / (and lower s slightly), suggesting that as the mechanical 
feedback enhances (larger a due to higher P c ), a smaller kick is able to trigger the flowing 
instability when the motors are susceptible. 

To verify the predicted dynamic phases and visualize the structural development, we per- 
formed dynamic Monte Carlo^ simulations on the model cytoskeleton as a "cat's cradle"— 1 ^. 
In these simulations we generated initially a three-dimensional random lattice of volume- 
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less nodes (mimicking the crosslinking proteins) and connected the nearest-neighbor nodes 
(defined by the first shell of the pair distribution function) with nonlinear elastic bonds^ 
(mimicking the filamentous proteins) at a given probability P c which characterizes network 
connectivity. Thermal steps obey Brownian dynamics^ whereas chemical moves follow the 
stochastic process defined by the model kicking statistics (equation H]). We postpone for a 
more detailed article describing all the simulations in different parameter regimes and only 
illustrate some qualitative features hereafter. The simulations reveal an interesting inter- 
play of network connectivity with the motor susceptibility dramatically affecting structural 
development. 

At a relatively high network connectivity, P c ~ 0.5 (average coordination number z ~ 6), 
force transmission through the bonds is efficient and the network structure remains statis- 
tically homogeneous in the presence of the motor-driven processes. Nevertheless, varying 
motor susceptibility drastically changes the dynamics. 

fluidized state: Under completely adamant kicks (s = 0) with a moderately large step 
size (I > /*), nodes rapidly become fluidized. Elastic stretching of the bonds imposes no 
constraint on the node motion resulting in vanishing localization strength and zero net flow 
(m = 0). Consequently, as shown in Fig. [2^, almost all the initially floppy bonds (in green) 
get stretched (in red) and the network becomes very tense. 

flowing state: At the other extreme, however, under susceptible kicks (s = 1) with above- 
threshold step size I > / t h, a spontaneous and self-sustained flow develops and the nodes 
vibrate about a steadily moving fiducial lattice (see Supplementary Fig. S2 for statistical 
characteristics of the flowing state). We found spontaneous flow both for regular lattices and 
for random structures. A plausible physical mechanism is as follows: disorder in the struc- 
ture, inherent in the quenched connectivity or dynamically generated through initial random 
motions, gives rise to local force asymmetry and thus local directed motion. Sufficiently large 
kicks then enhance an existing local force asymmetry and trigger dynamic instability of the 
quiescent state; the resultant nucleation and propagation of local coordinated motion, me- 
diated by force transmission and orchestrated by susceptible motor kicking, finally leads 
to a global concerted movement of the whole lattice. High motor susceptibility promotes 
cooperativity, resulting in "rigidity" of the structure and coherence of the collective motion. 
As can be seen in Fig. [2b, even though the system undergoes a flow instability, the network 
structure remains homogenous without significant local distortions, reflecting the enhanced 
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FIG. 2: Network structure for various dynamic phases, (a) Fluidized phase with s = 0. An 
initially relaxed (left) network rapidly tenses up (right) under completely adamant motor kicking. 
Localization strength of the nodes vanishes and there is no net flow, (b) Diffusive phase with 
s = 0.2. Spots of concentrated tense/floppy filaments are visible. Nodes exhibit enhanced diffusive 
motion with a finite localization strength. No spontaneous flow occurs, (c) Flowing phase with 
s = 1. Network structure remains homogeneous despite the spontaneous flowing motion, reflecting 
the enhanced rigidity of the structure and coherence of motion by susceptible motor kicking. L e = 
1.2, /?7 = 5, P c = 0.5 and I = 0.25. Red lines stand for tense filaments and green lines for floppy 
filaments. The small symbol at the bottom of each panel indicates the coordinate frame. 

mechanical stability and coherence of motion due to susceptible motor kicking. 

diffusive state: When the motors are only moderately susceptible to force (s : 0.2 — 0.5) 
yet are not sufficiently cooperative to drive spontaneous flow, the system exhibits enhanced 
diffusive relaxation toward the effective equilibrium characterized by T e g, leading to a homo- 
geneous network structure with modest local density fluctuations. The magnitude of density 
fluctuations and the tenseness of the network depend on the susceptibility. At relatively low 
susceptibility (s < 0.3), homogeneously distributed spots of concentrated tense or floppy 
filaments are visible (Fig. Wp)] as the susceptibility rises (s ~ 0.5), density fluctuations 
get weaker and the network becomes more homogeneous with a lower degree of stretching, 
closely resembling the flowing state (Fig. Therefore, for a homogeneous network at a 
relatively high connectivity the simulations verify the possible dynamic phases predicted by 
the analytical theory, that is for a sufficiently large kick step size, the system may undergo 
fluidized, diffusive and flowing states as the motor susceptibility increases. 

Another noteworthy dynamical and structural development occurs at intermediate con- 
nectivity P c ~ 0.3 (z ~ 3-4). In this situation network connectivity is sufficient for tension 
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percolation yet local force asymmetry becomes significant and widespread over the network. 
Under susceptible motor kicks with a considerable step size, dramatic spatial heterogeneity 
emerges and evolves into collapses and oscillations of the network in a particular spon- 
taneously chosen spatial direction (Fig. [3] and Supplementary Movie SMI). Apparently 
the overall tenseness of the structure is reduced by collapsing the network into clumps at 
the cost of a few highly stretched inter-clump filaments. Figures [3fc,d display the planar 
clumps where the floppy filaments (in green) concentrate and which are connected by highly 
stretched inter-clump bonds (in red). Collapse can occur in various directions and the fea- 
tures of the spatial patterns are independent of the system size. Our previous analytical 
mean-field study of an equilibrium nonlinear-elastic network— already suggests the possibil- 
ity of phase separation in this system at a finite effective temperature; the pressure exhibits 
a non-monotonic dependence on the node concentration leading to mechanical instability 
of homogeneous states. (Phase separation induced by contractile instability has also been 
predicted for active polar gels^.) When confined by boundaries, such as the cell membrane, 
these oscillating clumps may become stationary wave patterns with a characteristic length 
scale of modulation, reminiscent of the mitotic spindles. 



FIG. 3: Temporal development of the phase separation. At intermediate connectivity P c = 0.3 
and under susceptible motor kicking (s = 1) with a considerable step size I = 0.25, the initially 
homogeneous system (a) phase separates into oscillating clumps and voids (b-d). Notice that the 
relaxed filaments (green) become concentrated within the planar clumps as well as the presence of 
highly stretched inter-clump filaments (red). L e = 1.2, ^7 = 5. 

In the presence of susceptible motor kicking (s > 0.6), failure of force percolation at lower 
values of the network connectivity P c < 0.2 (z < 3) also yields phase separation but without 
a preferred direction of motion (Supplementary Movie SM2). Conversely when there are 
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too many bond constraints at P c > 0.6 (z > 7) there are significant mechanical barriers for 
motor dynamics which seems to slow down flow initiation and to reduce flow velocity. 




FIG. 4: Dependence of growth rate upon localization strength at different motor susceptibility. 
The 2D surface and the contour map of the growth rate r^/k 2 for flowing instability are displayed 
over the parameter plane extended by the localization strength a of individual particles and the 
motor kick step size I. Dq = 0.1 and k = 20. (a) s = 0.5; (b) s = 0.55; (c) s = 0.6; (d) s = 0.7. 
For susceptible motors with s > 0.6, the growth rate develops a non-monotonic dependence on the 
localization strength at a sufficiently large kick size, suggesting an optimal strength of mechanical 
feedback for an efficient flowing motion. 

The analytic stability analysis leads to a similar connection between the growth rate of 
the dynamic instability to the number of bond constraints via the localization strength of 
individual nodes. In Fig. 2] we display the two dimensional surface as well as the contour 
map of the growth rate r^jk 2 given by Eq. ( ITU]) in the flowing regime (r^ > 0) as a function 
of localization strength a and kick size I for a series of motor susceptibility s. Close to 
detailed balance, i.e., s = 0.5 (panel a), flowing instability emerges only at very high a 
and large /, and the growth rate increases with a. At s = 0.55 (b), a plateau in the 
growth rate develops at relatively high a range. For susceptible motors with s > 0.6 (c 
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and d), the growth rate exhibits a non-monotonic dependence on a, indicating that as the 
localization strength of individual constituents increases, growth of the flowing instability 
first speeds up and then slows down; in other words, there exists an optimal localization 
strength (or network connectivity) for the most efficient flow. This is consistent with the 
numerical observation that in a tense network (high a) that exhibits spontaneous flowing 
motion driven by susceptible motor kicks (s > 0.6), as the network becomes increasingly 
connected, flow initiation as well as the flow velocity first increases (for P c : 0.4 — > 0.6 or 
2:5—7-7) but then decreases (for P c : 0.6 — > 1 or z : 7 — > 12); nevertheless a flowing 
instability persists under these conditions (Supplementary Movie SM3-7). 

In principle motors could have slip bonds^i^ such that an applied force shortens bond 
lifetimes by lowering the energy barrier rather than slowing up-hill moves. This slip-bond 
behavior leads to a negative motor susceptibility which in turn leads to a negative effective 
temperature. This implies an intrinsic thermodynamic instability. We investigated the 
course of structural development in this thermodynamically unusual situation. We consider 
the case for s u = — l,Sd = where motors are insensitive to energetically downhill slope 
while they run faster when they go up against obstacles. Starting with a disordered structure 
at a relatively high connectivity (P c = 0.5), the motorized network rapidly develops into a 
highly ordered and tense structure as shown in Fig. [5l 



FIG. 5: Aster-like patterns/Bundle-connected poles formed under kicks of uphill-prone motors 
with a negative susceptibility. s u = —1, Sd = 0; L e = 1.2, /?7 = 5, P c = 0.5, and I = 0.25. (a) and 
(b) are snapshots of the system taken at the same instant from different view angles. 

This interesting behavior is not hard to understand: due to the propensity to capture 
uphill slopes (negative s u ), consistent with the negative effective temperature, the motor 
kicks try to maximize the total energy by separating the bonded nodes as far as possible 
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from each other; within the periodic cubic box the nodes are therefore concentrated at the 
corners giving the maximal inter-node separation and thus the highest degree of stretching. 
The resulting "aster-like" patterns closely resemble those formed by in vitro reconstituted 
active gels 3 - 3 -, where unidirectional movement of myosin II motors along the polar filament 
track toward the aster core (concentrated 'plus' ends of actin filaments) results in consider- 
able stress accumulation at the center, giving rise to the so-called "novas of asters". This 
is another demonstration of how motor susceptibility dramatically affects development of 
nonequilibrium structures. 

In sum we have derived an analytical expression for the stability limits of quiescent active 
gels and proposed a mechanism for spontaneous collective motion within a unified theoretical 
framework, based on a systematic expansion of the motor-driven nonequilibrium dynamics 
in terms of the kick step size. Simulations of a model cytoskeletal network further highlight 
that the interplay of network connectivity with motor susceptibility dramatically affects the 
formation of nonequilibrium structures: force percolation and mechano-chemical coupling 
conspire to drive and maintain the spontaneous flow, whereas adamant motor kicks generally 
promote fluidization. Significant force imbalance sensed by susceptible motors induces phase 
separation into oscillating clumps. Uphill-prone motors with a negative susceptibility give 
rise to a system at a negative effective temperature. Such motors drive the formation of 
aster-like patterns resembling what is seen in reconstituted active gels. 

We gratefully acknowledge the financial support from the Center for Theoretical Biolog- 
ical Physics sponsored by the NSF via Grant PHY-0822283 and a critical reading of the 
manuscript by Olga K. Dudko. 
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Supplementary figures 




Fig. SI Dependence of threshold kick step size upon mechanical gradient a and 
motor susceptibility s. We plot f(l,s;a,a) given by equation (11) in the main text 
as a function of kick step size I for k — 20, D = 0.1. The crossing point of /-curve with 

1(3 



/-axis locates l t h- (a) From right to left, a = 0.1,1,10 and 100; s = 1. / t h decreases with 
increasing a value (dashed arrow), indicating that strong mechanical feedback facilitates 
instability onset, (b) For a moderate a value (a = 5), no instability onset occurs for 
s < 0.5; only after s reaches approximately 0.6 does an instability emerge and increasing s 
(dashed arrow) leads to lower 1^. Meanwhile, for a given I, f increases with s, suggesting 
that more susceptible motors drive faster flows. 
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Fig. S2 Statistical characteristics for the flowing state. We plot the ensemble- 
averaged temporal correlation of displacement fluctuations defined via the correlation func- 
tion C S f{t COII ) = d,T[j^ J2f=l ^i( T )][jf <^( T + 4 corr)], where 8fi(t) = fi(t) ~ 

(fi(t)) with the bracket indicating an average over a steady-state time window. The integral 
averages the starting time r over a wide time range (to, tf). Correlations in three orthog- 
onal spatial directions (x: red; y: green; z: blue) all exhibit periodic oscillations in time 
(in Monte Carlo units) reflecting the spontaneous and coherent flowing motion across the 
periodic simulation box. Other correlation measures exhibit consistent oscillatory patterns. 
L e = 1.2, /9 7 = 5, P c = 0.5; / = 0.25, s = 1. 
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